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We propose a scheme to reconstruct the state of a two- 
mode Bose-Einstein condensate, with a given total number 
of atoms, using an atom interferometer that requires beam 
splitter, phase shift and non-ideal atom counting operations. 
The density matrix in the number-state basis can be com- 
puted directly from the probabilities of different counts for 
various phase shifts between the original modes, unless the 
beamsplitter is exactly balanced. Simulated noisy data from 
a two-mode coherent state is produced and the state is recon- 
structed, for 49 atoms. The error can be estimated from the 
singular values of the transformation matrix between state 
and probability data. 
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In quantum physics, a system can be completely de- 
scribed by its quantum state. On the other hand, in 
a given experiment one can only measure the probabil- 
ity of a quantity, and the process of doing so usually 
destroys any knowledge of the complementary aspect of 
that quantity. Thus to reveal the whole state a variety 
of measurements must be made on an ensemble of iden- 
tically prepared systems. Even when the systems are not 
identical one can recover a mixture of quantum states 
in the form of the density matrix. Optical homodyne 
tomography (OHT) is one such method of revealing the 
density matrix, tailored to determine the state of an elec- 
tromagnetic field mode jjj . In that case the availability of 
coherent states from lasers allows the measurement of the 
probability of a quadrature value, for a given phase differ- 
ence between the coherent state and the signal state. It 
also lends itself to a simple geometric interpretation: the 
measured distributions are the "shadows" of rotations of 
the Wigner function. Recent work has shown that it may 
be easier numerically to reconstruct the density matrix 
in the number-state basis directly from the distributions, 
without first going through the Wigner function M. 

In the field of atom optics, breakthroughs in the evap- 
orative cooling of atomic vapors have resulted in obser- 
vation of Bose-Einstein condensates (BEC) We may 
view these as the analog of a single-mode optical field, 
except that the field now has a finite mass and a quar- 
tic interaction term. The "semiclassical" aspects of BEC 
such as collective excitations and condensate shape are 
now a subject of thorough experimental investigation [Q. 
However, perhaps the most interesting aspects are due 
to the quantum states themselves. We might expect un- 
usual quantum states for several reasons. One is, the col- 
lisions are analogous to nonlinear optical susceptibilities 



which are known to produce squeezing |3|,|6|]. Second, an 
entangled state with a random relative phase can emerge 
between two condensates, when they produce spatial in- 
terference fringes in the atom counts |7j]. Also, in a closed 
system, quantum superpositions of different total num- 
bers are not permitted (although mixtures of different 
numbers are). For this reason, the OHT method cannot 
be applied to these BEC since coherent states are such 
superpositions of number states. A recent experiment || 
has demonstrated coherence within a single output pulse 
of atoms; however, the phase varies from one pulse to the 
next, and the exact nature of the output state is not yet 
known. Although theoretical work on atom lasers sug- 
gests that coherent atom fields will eventually be avail- 
able J! , we instead take a different tack and reconstruct 
the density matrix from a two-mode system with fixed 
atom number. 

We describe an atom interferometer, and a set of mea- 
surements sufficient to reconstruct the state. Using an 
angular momentum formalism suitable for the two-mode 
states, the reconstruction process amounts to the solu- 
tion of a linear algebra problem. This is accomplished 
with standard Fourier transform and least-squares meth- 
ods. Finally we reconstruct a two-mode coherent state 
from simulated noisy data, and compare the actual error 
with an analytical estimate. 

We consider a fixed number of atoms which have been 
cooled to below the critical temperature for BEC. The 
atoms are prepared in the two-mode quantum state 
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that is to be measured. The modes may correspond to 
e. g. different hyperfine states; such a BEC has just 
been produced in an experiment |[o]]. One of the modes 
has its phase shifted relative to the other by (f>, that can 
be varied. Then an operation corresponding to that of 
a beam splitter is applied to the two modes. The num- 
ber of atoms in one of the output modes is counted with 
the detector (see Fig. |l[). If the initial state preparation 
is considered to be a generalization of a beam splitter 
(which entangles a product state), the interferometer re- 
sembles a Mach-Zehnder arrangement. To collect data 
for tomography, for each of K phase shifts, we repeat the 



1 



experiment r times, resulting in probability distributions 
of counts at each value of 4>. 

At this stage we are making several idealisations. We 
neglect the collisions between the atoms, especially at the 
beam splitter where this would result in mixing the two 
modes. The beamsplitter is assumed to be lossless, with 
transmission cos 2 9. The action of the phase shifter and 
beamsplitter together on the density matrix p for the two 
modes is given by 



Pom = u{8, 4>ypu(6 



with 
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Finally we assume the detector has unit efficiency. This 
assumption is not crucial, as we can accumulate the same 
set of data when the efficiency is less than unity by count- 
ing atoms in both modes and keeping only the data for 
which the sum of counts was N. 

One possible realisation of the interferometer is an ex- 
tension of the recent output coupler built by the MIT 
group et al || . After the condensate is prepared in a su- 
perposition of \gmp) and \gvnp + 2) hyperfine states, the 
atoms are dropped out of the trap. The operation U{9, </>) 
is accomplished via two off-resonance optical pulses, con- 
necting the states in a Raman transition, which also in- 
troduces a relative phase shift due to the phase between 
the pulses. Finally an on-resonance rf pulse transfers 
atoms from the untrapped \gmp) state to the trapped 
\g'm,F) state; the separate groups of atoms can then be 
counted. 

To simplify the notation and the calculation of matrix 
elements, we use the formal equivalence between the al- 
gebra for two harmonic oscillators and that for angular 
momentum |nj . We write the state 

\n)i (g) |JV - n) 2 = \j +m)j. ® Ij — m) 2 as \m) (5) 

where j — N/2 and m — n — j . The 2j + 1 states \m) 
have all the properties of the eigenstates of J 2 and J z 
and important operators are 
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2«1 



J- = ala 



ai — a\a.2 



(6) 
(7) 

(8) 



The combined effect of the beamsplitter and phase shift is 
seen to be a rotation by —29 about an axis = x cos <f>— 
ysin<^, 

U(9, 4>) = exp [t0(J+e** + J-e**)] 

= exp [i29J ■ n ] . (9) 

An explicit expression for the number-basis matrix ele- 
ments of the rotation is 
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where the sum is over all values for which the arguments 
of the factorials are nonnegative. 

The probability of m counts at the detector, for a phase 
shift setting of <j> is 



P m (<j>) = {m\U(6,4>y Pin U(6,<p)\m) 
or in the number-state basis 
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This last equation shows that T is simply a linear opera- 
tor taking the density matrix (an (N + l) 2 element com- 
plex vector) to the probability data (an (N+1)K element 
real vector); we wish to invert this operator to recon- 
struct the density matrix. This is formally done by defin- 
ing the positive semidefinite, hermitian (7V+l) 2 x (iV+1) 2 
matrix 



M = T f T, 
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where means the conjugate transpose of the matrix 
T, and computing 
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provided M is nonsingular. This is actually the best 
solution in the least squares sense Q . 

For the value 9 = j, corresponding to a balanced 
beamsplitter, the matrix M is singular. Since this might 
be somewhat surprising we prove this result, and give an 
example of two states which will produce identical prob- 
ability data. First we note that if there exists a vector 
v such that Mv = 0, this implies Tv = 0. Making the 
Ansatz 



v = f(J z ) 
/(-<?) = -/(<?), 
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and using the formula for T in the number state basis 
Eq. (H|) we find 
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Then, using the fact from Eq. ((TTj) that 
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yields zero for the sum of an odd function in Eq. (18), 
for any values of I and <j>. This proves that the matrix 
M is singular for 9 = j. To illustrate this physically, we 
choose the pure states 
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having the respective density matrices 

Pi = |"i>(m| 

P2 = | — m)(-m|. 
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Then the difference vector of the density matrices is in 
the nullspace of T: 

{1 if x = m 
-1 if x = -m (24) 
otherwise 

so that Tpi = Tp 2 - Clearly the balanced scheme cannot 
distinguish these two states because there is an ambiguity 
about which number of atoms belongs to which mode. 
This difficulty does not arise in OHT because in that case 
one of the inputs is already known to be a coherent state. 
Of course, varying 9 will enable these two situations to be 
distinguished, but we would like to keep the experimental 
set up the same for all of the data collection, so we require 

The problem of solving Eq. ([12]) is considerably simpli- 
fied by taking the Fourier series of the probability data, 
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(In practice, this is approximated by the discrete Fourier 
transform on a set of data for many <fi equally distributed 
about the circle.) By taking succesive values of r from 
to 2j, each diagonal of the density matrix can be solved 
for independently. Since all the operations are linear the 
calculations are conveniently done in MATLAB. 

We will see below that for elements far from the main 
diagonal, the data can be insufficient to accurately re- 
construct them. Since the trace of p 2 is less than unity, 
each norm of a diagonal is also less than one. Taking this 
prior information into account, we can reduce the norm 
of each diagonal by using Tikhonov regularization jl3j . 

We present a simulation of the probability data and 
reconstruct the density matrix from it, for the pure state 
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This state is the projection of coherent states for modes 1 
and 2 onto the subspace with fixed total number A" ; it 
can be produced by an rf pulse as in || . The parameters 
chosen were N — 49, $ = 0.54, (p = 0.13. To simulate 
the noise associated with calculating probabilities from 
histograms, we add to each probability a term 



SP m (<l>) 
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where r is the number of trials for that value of phase 
and g is a Gaussian distribution with zero mean and unit 
variance. In the example below we take r = 2000 and 
K = 180, i. e. 2000 trials every 2 degrees. In Fig. | 
we plot probability of m counts at each phase shift. The 
original and reconstructed density matrices are compared 
in Fig. ^. The reconstruction is most accurate (in both 
magnitude and phase) for the elements close to the main 
diagonal. 

The error in the solution to the least squares problem 
can be found using the singular value decomposition of 
each matrix T |l2|] . Specifically, if the norm of the error in 
F(r) is ||<5F(r)|| then the norm of the difference between 
the original and reconstructed r-th diagonal is 



l«Pr|| = 



\ 



E 



Pi,i-\-r 



2 < 



\SF(r)\\ 



(29) 



where a is an eigenvalue of yT^T. Using the definition 
Eq. (25) to transform the variances SP m , and the fact 
that probabilities for a given phase must sum to one, we 
find 
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Thus our upper bound is independent of the actual data, 
depending only on N, r and 9. In Fig. ^| we plot 
the actual and estimated norm of the error, divided by 
\/N — r + 1 so that it is per element, as a function of 
the distance r from the main diagonal. We have not 
been able to find a simple formula for the eigenvalues 
a = cr(N, r, 9), but the general trend is for smaller mini- 
mum values as r approaches N. The most difficult states 
to reconstruct will be highly entangled ones, which de- 
pend on the high-frequency components of the proba- 
bility. Provided one has some idea of which diagonals 
are needed this analysis is useful in choosing the optimal 
transmission. 

A different type of error might occur if the number 
of atoms originally trapped is not known precisely, or 
(equivalently) if losses occur between the preparation of 
the original state and the detection of atoms. We are un- 
able to quantify such errors, since states of fixed number 
are assumed. For mixed states with narrowly distributed 
total numbesr, but no fine-scale oscillations in density 
matrix elements, we can expect the above algorithm will 
produce a qualitatively correct picture. 
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To conclude, we have proposed an atom interferome- 
ter for measuring the quantum state of a two- mode BEC. 
An algorithm for calculating the density matrix elements 
was applied to simulated probability data and success- 
fully reconstructed a sample state. The errors can be 
estimated before data collection, and may be minimized 
by choosing an appropriate beamsplitter transmission if 
some qualitative features about the state are known. 
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of the Royal Society of New Zealand and The University 
of Auckland Research Fund. Frequent discussions with 
Tony Wong and Michael Jack provided many insights. 
We would like to thank M.-O. Mewes for his detailed 
suggestions regarding the proposed interferometer. 



FIG. 1. Schematic of atom interferometer needed for state 
reconstruction. One of the modes acquires a phase shift <j> 
before they both undergo a beamsplitter operation (BS). An 
atom counter (D) is placed at one or both outputs. 

FIG. 2. Probability as a function of number of atom counts 
at one port and relative phase shift between modes. Com- 
puted from Eq. (fi^) with simulated noise Eq. (|28|). 

FIG. 3. Density matrix of (a) the original state and (b) the 
reconstructed state, in the number basis. Magnitude is height 
and phase is color. 
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FIG. 4. Data points: Norm of the difference, per element, 
in the original and reconstructed density matrices, as a func- 
tion of distance from the main diagonal. Solid curve: Es- 
timated norm of the difference, per element, based on the 
singular values of the transformation matrices. 
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Fig. 1: Atom Interferometer 



Bolda, Tan & Walls Fig. 2: Probability data 




Bolda, Tan & Walls Fig. 3(a): Original density matrix 




Bolda, Tan & Walls Fig. 3(b): Reconstructed density matrix 




